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Abstract 

System reliability estimates are generally made using a mode}! 
which assumes independence between components, and results are often 
claimed to be conservative. For a two component serial system bivar. 
iate distributions are developed for three cases: (1) bivariate expon 
ential, (2) bivariate geometric, and (3) a composite exponential, 
geometric bivariate. These distributions are then utilized to inves- 
tigate the reliability of a two component serial system when an 
estimate of the correlation coefficient is available. The estimate 
of the reliability thus obtained is then compared with the corresponding 
estimate obtained by use of the model which assumes that the system 
reliability is the product of the component reliabilities. The dif- 
ference between these two estimates is tabulated for values of the 
correlation coefficient between -.25 and +.25, for each of thet" -e 
bivariate distributions. Conditions under which the effect is maximum 
are explored and a method of approximating the reliability difference 


is suggested. 


The writers wish to express their appreciation for the assistance 
and encouragement given them by Professor Walter Max Woods of the 


U. S. Naval Postgraduate School in this investigation. 
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Section l 


INTRODUCTION 


It is known that various types of interdependence can exist 
between components of a general system. The resulting effects on 
the system reliability will be a function of the correlation between 
the components considered. A practice in reliability studies has 
been to assume that ignoring correlation effects would lead to a 
conservative estimate of the reliability. We propose here to study 
the validity of this assumption for various values of correlation 


and further to study quantitatively the resulting reliability estimates. 


To narrow the scope of the problem, we have considered a serial 
System of two components and examined three distributions; (1) bivariate 
exponential; (2) bivariate geometric; and (3) a composite bivariate 


distribution where the marginals are exponential and geometric. 


Each of these distributions is examined separately. A summary 
of their characteristics is given in Table 1.1. A comparison is made 


of the reliabilities defined by 


li 


Beeic eae De Sate oS tl = ke Ler ay = Ge 


pix2t, yve2tJ (122) 


Ry (t) 


for positive and negative values of 0. 


The resulting estimates of reliability are analyzed and the 
effects of correlation on system reliability are evaluated. Exact 
confidence limits are derived where possible and approximations are 


considered otherwise. 





Table 1.1 


Probability Distributions 


Exponential Family 


Density : f_(t) = ie exp (-t/a) , t 20 
le a 
0 t= 0 
Parameter ; a 
a > 0Q 
Mean : a 
Variance : a 
Geometric Family 
Mass function; p,(k) * os (l-p) ele 58 
K : 
0 otherwise 
Parameter : Pp 0 Spsil 
Mean : p/l-p 
Variance : ay Clee) 
Bivariate Exponential 
Density : fy yy) = £, (x) f(y) li +v (2F, (x)-1) (2F,,(y}-1) J 
where f(x) = ep (=x/a) 
X a 
1 x — 0 
f(y) = 5 exp (~y/b) ee 
-] Sv s1 
FY (x) = 1 - exp (-x/a) a20Q 
b20Q 


Fy(y) = 1 - exp (-y/b) 


Parameters : a. bev 





Section 2 


GENERAL CONCEPTS OF RELIABILITY 


2.1 Definition 

Reliability, as the term is used in mathematical statistics, 
has exact meanings. It can be calculated, objectively evaluated, 
tested, and designed into equipment. A definition, as given by Lloyd 
and Lipow les is "the probability of a successful operation of the 
device in the manner and under the conditions of intended use". 
Mathematically then, reliability, at some point Xo» is the probability 
that a random variable, X, representing the operating life of some 
device, will equal or exceed the point Xo: It can be seen that if 
the probability distribution of the random variable is known, then 


the reliability at any point may be calculated. 


A function representing the operating life may be either a con- 
tinuous or discrete type of probability distribution. The interval 
of operation can be thought of as a time interval, in which case the 
probability distribution will be continuous, or as a number of oper= 
ations (such as the turning on of a switch or relay) in which case 
the distribution is of the discrete type. Thus, for a continuous 
distribution the random variable T, to be considered, is the time 
to failure, while for a discrete distribution the random variable might 
be the number of power turn-ons (PTO's), K, before failure. Combin- 
ations of these are also possible. For example, the number of starts 


1. Square brackets refer to correspondingly numbered references 
Shown in Bibliography. 





of a jet engine and the duration of operation may determine the lite 


of the engine. 


In the continuous case, the reliability, R(t), of a component 
is the probability that the component will operate at least for some 
tame, tt... Thus, if f(t) is the probability density function for T, 
then reliability or probability that the time to failure will exceed 


or equal t, is 


te @) 


R(t) = P[tT2t] = feqceree (22) 
C 


In the discrete case the life of the component may be measured 
by the number of PTO's prior to the first failure, where the random 
variable, say K, may take only discrete values 1, 2, 3, ... . To 
determine the reliability at some fixed value k, where R(k) = Prob- 
ability [the number of PTO's prior to first failure will exceed kJ, 
we must determine the probability that K will take any value greater 


than k. 


If the random variable K has a certain probability mass function 
@ 
p(k) such that 0 $ p(k) $1 and ay p(k) = 1 where p(k) = P Lk = kd, 
and the probability, p, that a certain single trial results in a 
success, is constant for any particular trial, ie, each PTO is a 
Bernoulli trial with fixed probability,p, of success, then the mass 
function is p(k) = ee - p) (the probability of k successes and 


then a failure). 


The reliability at some point Ke is then the probability that 
the random variable K will take any value greater than or equai to kK, 


Ores 





Hence: 


8 


Rae = = p”q (2.2) 
where q®=il-=-p 


Throughout this investigation we will assume that each triai 


does constitute a Bernoulli trial with fixed parameter p. 


2.2 Component Interdependence 

In determining the reliability of a system consisting of components, 
two basic configurations are of interest (1) a series system and (2) 
a parallel system. The series system consists of components put to- 
gether in such a way that each component must operate for the system 
to operate, while in a parallel system the system will operate if any 
one component functions. This investigation will consider only a 
series system consisting of two components, with the probability 
distribution of each being known and the parameters of the distri- 


butions estimable for each. 


If it is assumed that there is no component interaction in a 
series system, then the distributions are statistically independent 
and reliability of the system can be determined from the reliability 
of the components by multiplication of the component reliabilities. 
It should be noted, however, that the product of the component reli- 
abilities may, under certain circumstances, indeed yield the system 
reliability even though strict independence does not hold. If there 


is interaction between components, then the product rule does not in 





general hold. Rosenblatt [16] states that frequently the observation 
is made that "assessments of system reliability based on the assumption 
that component failures occur independently of one another are approx- 


imate and usually excessively conservative". 


We will show that interaction between components can reduce as 
well as increase the reliability of a system. An example of inter- 
action which reduces system reliability could be the case where two 
components, each of which produces heat when operated, causes a 
temperature environment which reduces the life of one or both com- 
ponents and hence of the system. Each component when operated 
separately, however, might produce less heat than is required to 
affect the life of that particular component. On _ the other hand in 
certain electronic devices a high temperature environment might be 
beneficial, in which case the interaction described would enhance the 


reliability of the system. 


2.3 Statistical Estimation 

In the estimation of reliability of a system from information 
available on the components, if the probability distribution family 
is known then the problem reduces to determining or estimating values 
of the parameters of the distribution. Two methods are available for 
doing this:(1) point estimates and (2) confidence interval estimates. 
A point estimate is the value of a statistic based om some experi- 
mental measurements. An example is the maximum likelihood estimate 


(M. L. E.), which has certain optimal qualities, see Mood Liga 





A two-sided confidence interval estimate of a parameter is 
obtained by selecting two random values L and U such that, given 4 
number 0 < yY < 1, the statement that the random interval om ul 
covers the parameter may be made with probability 1 - Y. A one- 
sided confidence interval is obtained by selecting a single random 
value L where the corresponding confidence statement is, the interval 


[L, ©] covers the parameter with probability 1 - Y. 


The method by which experimental data is obtained is called the 
Sampling plan. Several sampling pians are available for estimating 
parameters of the exponential distribution. Those considered in 
this investigation are: 

Sampling Plan I: test N items of a given type until all fail 


and observe the N times to failure. 


Sampling Plan II: test N items of a given type until r of them 
fail, where r SN is fixed before starting 
the test. The observations are then the first 


r failure times. 


Sampling Plan III: test n items of a given type to a preassigned 
time to: Let r denote the random number of 
failures in the specified time and observe 


the r failure times. 


The only sampling plan which will be considered for the geometric 
case is as follows: observe the number r of failures for a pre-assigned 


number, N, of PTO's and any number of items of the given type. N is 





fixed in advance of the test and r, the number of failures, is a 


random variable. 


Association between two variables is estimated by several co~ 
efficients, the most notable of which are; 

(i) Product moment correlation coefficient, 

(it) Spearman's rank correlation coefficient, 


(iii) Kendall's 7T correlation coefficient. 


The performance of these correlation coefficients for general 


bivariate distributions are compared by Farlie Els 


2.4 Summary 

In summary, our procedure for estimating reliability might be 
broken down into three steps: (1) to establish the type of statis- 
tical distribution which describes the failure phenomenon, here we 
assume this known, (2) to estimate the parameters which completely 
define the distribution, and (3) to utilize the knowledge of the 
distribution with the estimates of the parameters to estimate the 
reliability. In many present day applications a simplified model, 
the independent serial system model, in which the system reliability 
is calculated as the product of the reliability of the components, 
is used even though the actual reliability may be quite different 
and not always greater. Jt takes but little reflection to realize 
that in many applications an underestimate even by a small percentage 
might have great consequence on the cost of a large development 
program. It therefore seems very appropriate to attempt at least a 
start at a procedure by which the reliability can be estimated taking 


account of any component interaction. 


=Ja 





Section 3 


BIVARIATE EXPONENTIAL 


3.1 Derivation 

In this section we shall be concerned with developing a mathe~ 
matical structure to represent the time to failure distribution of 
a two component system where each is exponentially distributed. 
Using this structure we shall determine the reliability of the system 
with known correlation and compare this with the product rule system 
reliability. We shall consider the marginal density functions to be 
of the form 


| 1 | 
f(t) = — exp (-t/a), t=" 0. ost) 


The corresponding cumulative distrubution functions are of the form 


t 


F(t) = [ro dt = 1 - exp (-t/a). (3.2) 


Gumbel [3] has demonstrated a general method for deriving a 
bivariate distribution function and the corresponding density function 


from two known distribution functions by applying the formulas 
Fyy06y) = Fy) Fy) [1 + v G = FG) C1 - FyG))] (3.3) 


where - 1 Sv $1 and consequently, 


fyyOoy) = &00 fy) [1 +v FG) = 1 QFG) - D (3.4) 


He has further derived two specific bivariate exponential distribution 


functions each of which is restricted to ranges of correlation less 


ene es 





than unity. We have chosen the more symmetric of the two distri- 

butions and have used it exclusively. For this case, where FY (x) 

and FY Cy) are both exponential, the bivariate exponential distri- 

bution and density functions become 

Fyy(x,y)= [1-exp(—x/a ) l-exp(=—y/b ) [ 14vexp(—x/a-y/b )| 9 (3.5) 
Oey oe © 


and 


1 
fy, (xy)= ap exp (=x/a~y/b)| l+v | exp (~x/a) -1] 
[2exp(~y/b) ]; (3.6) 


These can be shown (see Appendix A.1) to possess all the required 


properties and in particular the correlation coefficient can be ex- 
pressed as 

p= v/4. Sys. (3.7) 
From this it is seen that the correlation of this particular distri- 


bution is restricted to the range ~.25 SP Ss .25. 


3.2 System Reliability 
As previously asserted, reliability is a probabilistic statement 
about the operating life of a unit. If we define reliability in symbols 


as 


R(t) = Pit2c] en 


a ie 





then for the exponential case we have 


co 


R(t) = [ire dt = exp (-t/a). (3.9 
E 


SZ 


If we now consider the reliability of a system of two expon-~ 
entially distributed serial components, we express the system 


reliability as 


co co 
R(t) = P [fx2t, Y2t]= a [ fy (x,y) dx dy. (3.10) 
t 


We seek to investigate the consequences ot the assumption that 
RCE) euR Ct) one Lx2t] ply2cl] (3.11) 
when something is known of the correlation, Pf. 
Using the product rule, eq. (3.11) is indeed true. However, 


for the general case the system reliability is found (see Appendix A.1) 


to be 


R(t) =exp(-t/act/»)| 1+v[2-exp(-t/a)-eno(-t/e)texn(-t/a-t/> ] : (3512) 


This is monotone increasing in — and hence is restricted to the values 
of P between -.25 and .25. To obtain a more complete analysis the 
system reliability was derived for P = 1 by assuming ome one variable 
was a linear function of the other, say Y = cX where c > 0. In this 


special case the system reliability is given as 
R(t) =P [ x=t, y=] = vax| P [x=t], P [x=*/<]] (3.13) 
since for c $1, R(t) = exp (-t/a) and for c > 1, R(t) = exp (-t/ca). 


als 





To obtain a quantitative expression for the effect of correlation 
on the system reliability,a reliability difference function was con- 
structed. To eliminate the need for considering both component 
parameters in this function we denote their ratic as s = b/a and 
use this as a single parameter. Also, to avoid considering specific 
values of the operating time t, we shall use the ratio of the oper- 
ating time to the mean time t/a as a parameter. With this notation 
we shall define the reliability difference function, denoted as AR(t), 
as the difference between the system reliability when there is cor- 
relation, R,(t), and the system reliability when there is no cor- 


relation, R(t). Hence 


A R(t )=Ro(t)-Ry (+ )=4/0 exp| -t(st#2 )/sal ji-wxr(-t/a )-exp(-t/sa )t+exp {-r(s#1)/sal 
(3.14) 


Equation (3.14) is valid for -.25 =p S .25, but for p= 1 it is 
AR(t) = exp (-t/sa) ee exp (-t/sa) J. (3.15) 


Values of the reliability difference functions defined in 
equations (3.14) and (3.15) are given in Tables 3.01 to 3.10 for 
values of the t/a = .05(.05)2.5 and s = .5(.5)5. Further, they are 
plotted in Figs. 3.1 and 3.2 with respect to t/a. The curves vary 
linearly with P up to .25 so only the curve for Pp = .25 is shown 
along with that for ? = 1. To obtain values of AR(t) for P other 
than .25 the tables carry these for P = .05 and P = .15. At other 


values of § a linear interpolation will provide the answer. 


« ]3- 





For s = 1, when the component parameters are equal, the curve 
attains a maximum at t/a = .69315 and it is here that the effect of 
correlation on the system reliability is a maximum. When the par- 
ameters a and b are unequal, the curve peaks at different values of 
t/a depending on the value of s. For P less than .25, the maximum 
value of AR(t) is appreciably less than one-fourth that for Pp = 1 
at any s other than one. The maximum value of AR(t) is .2500 for 
P = 1] and .0625 for p = .25 which is linear and consistent with 


previous knowledge. 


Example 3.1 


As an example let us compute the reliabilities and correlation 


effects for a system at time t = 500 hours where a = b 1000 hours. 
Here t/a = .5 and the component reliabilities are R(t) = exp (-.5) = 


-60653. The system reliability for the independent case is R(t) = 


.36788. From Table 3.02, AR(t) for t/a = .5 and s = 1 is .05695 


for Pp ® .25 and AR(t) # .23865 for P 21.0. For Pp = + 1.0 the 
product rule underestimates the actual system reliability by 65%. 
Even for Pp = + .25 the underestimate is 15.5%. 
3.3 Confidence Interval 

A comprehensive discussion of the concepts and procedures in 
deriving estimates of the parameters and confidence intervals for 
the resulting reliability estimate is given in Chapters 7,8, and 10 
of Lloyd and Lipow [iJ]. We shall consider three sampling plans which 


are most likely to be utilized in obtaining the failure data: 


(1) testing N items until all fail, (2) testing N items until r <N 





fail, and (3) testing N items for time tO and observing r, the number 


of failures. In each case the underlying distribution is exponential. 


3.3.1 Testing N items independently until all faii 

In this case, a sample of N items are put on test and the test 
is concluded when all have failed. The times to failure ty» paige t 
each measured from the time the item was "turned on", are recorded. 


The intended life of the item is t, and it is required to demonstrate 


the reliability with confidence level ( 1 = Y). 


A maximum likelihood estimate for the parameter of the exponential 


distribution is 


oe | 
a = 7 a ts x (3.16) 
1=] 


It can be shown 11] that the function C. = 2N @/a has the chi-square 
distribution with 2N degrees of freedom. Hence a lower confidence 


limit on R may be derived from the expression 


pl2n&/a> X? ae Relay (3.17) 


This lower confidence limit, denoted by L , is 


L = exp (-t [ 2n 3 3.1 
pk X pee, N @) Fe LE 


For the case of two independent exponential distributions E08) 


and f(y) > if we define the ratio of the parameters as s = b/a then, 


from (3.12), we can denote the reliability as 


R(t) = exp (-t (s+1)/sa). (3.19) 





Since C. = 2 ny 8/a has the chi-square distribution with 2 ny degrees 
N\ 
of freedom and Cy = a b/b has the chi-square distribution with 
=2n @/at2n b/b. 
,b x y 


has the chi-square distri- 


2 ve degrees of freedom, let us define C. 
Then, by the reproductive property, oe b 
ae 
bution with 2 ny +72 ay degrees of freedom. If we substitute b = 5s a, 


we get 


a A 3 901 
oP: a (n, sa+ a b)/s a. (3 20° 


Using the same analysis as for the univariate case, we see that 


P | 2058 + ny b)/sa > 


eae (3.21) 


2(nmy)2 1 ~ 


Hence a lower confidence limit for R is 


L = exp =t (sct1) Xa(mgtt) -y 


with confidence level l = Y. 


a(n shtn b (3.22) 
ny 5°) 


3.3.2 Testing N items until r <N fail 
In this case a sample of N items are put on test and the test 
is concluded when some predetermined number r SN have failed. The 


times to failure t t are recorded where each ts is measured 


Ae a 
from the time the item was "turned on''. There arises the consideration 
of the immediate replacement of failed items. If replacement is used, 


the number on test is always N, whereas in the nonreplacement case, 


the number of items on test eventually drops to N=- r+]. 


LG 





It can be shown [1] that in either the replacement or non- 


replacement case the quantity (2 r a. y/@ has the chi-square dis 


r 
tribution with 2r degrees of freedom, where 2 y t. + (N-r) t.. 
i 1=1 
Hence the procedures and results derived in section 3.3.1 can be 


directly applied, and we observe that 


= 2 A a 
L = exp [tx co 1.9f** a0] (3.23) 


is the lower confidence limit on R(t) in the univariate case at 
confidence level 1 - Y. Similarly, for the case ot two independent 


exponential distributions a lower confidence limit on R(t) is 


L = exp| =t(stl e a. b (3.24) 
xp} =t{ x a rgtry) ; yf? By tty yon ) 


at confidence level 1 = Y. 


3.3.3 Testing N items until time ee 

In this case we test N items for a predetermined time to: We 
note the number of failures, r, in this time and record the times to 
failure ty shate Ce These t. are random variabies with density 


function 


Je : 
3 = — 2 os < < aS mz 
f(t) A Bex (-t/a) 0 t te (3: Ze 


where A = 1 - exp (-t /a). This is a density function since when 
integrated over the range of t the result is unity and F(t) is non- 


negative. The MLE of a is then given by 


Ast lee oe 
QA = ae by + (Nor ) ‘| (3.26) 


als 





hence an estimate of the reliability is 
R(t) = exp (-t/4). C3127) 


In order to obtain a confidence interval, we note that the 


random variable r has a binomial distribution with parameters N 


and p = 1 = exp (-t /a) A. Hence we see that 


L = (1 - @ (r) a (3.28) 


where @ (r) is the solution of I [n-r, rt1] = 8 and I (a,b) is 


1-a(r) 


the incomplete beta function tabulated by Karl Pearson. 


For the bivariate case the procedures are extended in the 


previous manner. If we denote the estimators as a and a, where 
Ux 


Bees = 2. ty + (Nex) ty 


and 
2 ry 
= y + (No 
rya t + ( Ly) t ee 


vos 


then we have 


R(t) = exp [-t (1/2 + 1/8.) ]. (3.29) 


Returning to our previous example, if we were to test 5 items 
of each component and we estimated the parameters as a = 975 and 


A 
b = 1050 then a 95% lower confidence limit on R(t) would be 


L = exp (-500 ( X* 20:95)/10125) = .53585 (3. 30) 


for the independent case. As a possible approach to a solution to 


are 





the dependent case we might try a procedure such as this: (i) obtain 
the independent confidence limit as above, (ii) by using one of the 
standard procedures derive an estimate for 9, (iii) find AR(t) at t/a 
using the value of R and (iv) add algebraically this AR(t) to Lab 

to obtain the final estimator. No claims are made as to the accuracy 
of this method, it is merely suggested as a possible means of ob= 


taining a bound on reliability in the face of known correlation. 


3.4 Approximating the Effect 
If we examine the reliability difference function defined in 
section 3.2 and use the notation of eq. (3.12), we may rewrite eq. (3.14) 


as 


AR(t)=40 exp(~t/act/»)| 1=exp(-t/a)=-exp(-t/b)+exp(=t/a=t/b) |: (3.31) 


By regrouping the terms slightly, we can put this into the form 


AR(t )=4, exp (t/a) l-exp (=t/a)] exp(=t/b) l=exp(—t/b) (3.32) 


and we see immediately that this is merely 4P times the product of 


the component reliabilities and their unreliabilities for -.25 Sp Ss .25. 


As a check if we let t/a = .5 and t/b = .25 and p = .25 then 


AR (t) -04111 which checks with the result in Table 3.04 for t/a = .5 


and s 2. This was to be expected since in this particular case the 


approximation is identical to the exact method. 
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3.5 Summary 

For a two component serial system the effect of correlation on 
system reliability is greatest at t/a = .69315 when the parameters 
are equal and the maximum AR(t) is .0625 for Pp = .25 and .2500 for 
P = 1.0. When the parameters are unequal, the point of maximum 
effect due to correlation varies directly as the ratio of the 
parameters. For the restricted bivariate distribution the maximum 
AR(t) is less than one fourth of that for the case P = 1] for values 
of s other than one, and this is likely due to the bivariate dis- 


tribution used. 


Lower confidence limits have been defined for the independent 
case under three sampling plans. For the dependent case an "ad hoc" 


procedure is suggested. 


3.6 Description of Graph Format 

At the end of this section and the succeeding sections are located 
the figures referred to in the text. These figures are presented in 
a standard format except for the axis scaling. This is indicated by 
the legend directly under the figure title. The following example 


illustrates the notation used. 


Example 3.2 
An axis scaling legend of X AXIS SCALE = 2.00 E + 02 
Y AXIS SCALE = 1.00 E - Q2 
2 


is to be read as "the X AXIS is marked off in units of 2.00 X 10° and 
the Y AXIS is marked off in units of 1.00 X ae A general legend 
would read "the X AXIS is marked off in units of A X 10°" and be 
denoted as X AXIS SCALE = AE + B. 
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Section 4 


BIVARIATE GEOMETRIC DISTRIBUTION 


4.1 Derivation 
In this section we shall develop a discrete bivariate dis- 
tribution based on a general class of joint distributions. We 
shall take as marginals the geometric distribution and then evolve 
a relation from which the reliability of a two component system may 
be determined. Using this relation we shall then examine the 
effect on the reliability resulting from various values of correlation 


between the two components. 


To investigate the effect of interdependence of two components 
when the probability of successful operation of a component follows 
a discrete distribution, a model applicable to ''power turn-ons" was 
used. In this model the number of power-turn-ons prior to the first 
failure is considered a random variable, X. It is assumed that the 
probability of a success (ie, a switch functions properly, a relay 
opens, etc.) is a constant, p. The probability of a failure, q, is 
then (l-p). The reliability at k PTO's R(k), is then defined as the 
probability of at least k successes before the first failure. The 
random variable X then follows a geometric distribution: Py (x) = pq 


where the cumulative distribution is given by: 


xceL 


4 
P(x) =P[x2x]= Y pqg=l-p (4,01) 


i=0 


(see Appendix A.2) 


= 33- 








SS oe 


The reliability of the component, ie, the probability that the random 


variable X will equal or exceed some constant, koe is given by 
as => = wz, = Ky ae ko | 
R(k,) = P[XEk,] = 1 = P[xsk,-1]=1-[1-p°]=p (4.2) 


In order to Study two components, each having a probability dis- 
tribution as defined above, a joint distribution, having identical 
but not necessarily independent, geometric distributions as marginals, 
was needed. Such a joint distribution is not unique but a class of 
such functions was determined by D. J. G. Farlie [2] in order to 
compare various correlation coefficients. (This investigation will 
consider only the product moment correlation coefficient, which is 


the more well known of the various types). 


The joint mass function: 


Pay (9) = pxCepy(y) | 1v[2P(e)pyo)-2 ]| Reve] (409) 
x= 0, 1, 25 eee 
y = OQ. i 25 seo 


was used. It can be shown (see Appendix A.2) that 


(1) Pyy sy) =0 


oO co 


(2 “y) = 
) By Ly Pa! y)=1 


<x 


(3) 2, Pyy(%sy) = py(y) and )) py (x,y) = py(x) 


ane 





Also it can be seen that when v * 0 the joint distribution reduces 
to the product of the marginal distributions. Therefore, this 
Pyy (>) does meet the conditions for a discrete bivariate dis-~ 


tribution function. 


To determine the functional relation between the constant v in 
Pyy (*>¥) and the product moment correlation coefficient, P, the 


defining equation for P was used. 


p= B [xv] - = [x] E [y] (4.4) 


Ox Oy _ 


Using Py for the parameter of the first component, ie, the 





component characterized by the random variable X, and P, for the 
parameter of the second component, the value of the correlation 


coefficient, P , was determined as; 


Vv VPL Po 


ciate (405) 
(1+p,) (14p,) 


(calculations are shown in Appendix A.2) 


From this relationship it can be seen that the values of ( vary 
from 0 when P, Or Py are zero to a maximum of i when P) = Py * i 
Now since -1 Sv $1, which implies that -% S$ P S 4, and hence this 
joint distribution function is satisfactory as a method by which to 
examine the effects on the reliability of the components having a 


correlation within this range. 





With this joint mass function, an equation for the reliability 
R(k) =plxe ko Ye kJ, was developed (see Appendix A.2). Using 
the previously developed relation between Pf and v, the reliability 
as a function of the correlation is: 


k, k, Pl1+p,)(14p,) k 
ce 


If it is assumed that the product rule holds, then the system 
reliability is the product of the component reliabilities. 
‘ k Kk. =*, k 
R, (ko) = Py} 9 Py oO ™ (PyP5) © (4.7) 


The difference between the two reliabilities, AR (Kk); is then 


the difference Ry (ko) -R, (ko): 


4.2 System Reliability 


To study the reliability difference the functional relation for 


AR (k.) was determined. 


- /A(i+pz) (145) 


PyPo 


AR(k 
: 2 2 


Ko 2k k 2k 
[P Oosawa © (58) 
eek 
Using eq. 4.8 and assuming P) ™ Po» the derivative with respect 
to ky of AR (Kk, ) was set equal to zero in order to determine the value 
of Ky at which the function attained a relative maximum (development 


shown in Appendix A.2). The value of ky at which this maximum occurs is 


aya 





in 2 





ae ee (4.9) 
For p = .999 this gives a K for maximum difference of 693. Taking 
the mean life of a component as E and taking the ratio of Ko to 
this mean life, gives k /P/4 = .693 ™ .7. Note that for the case 
of the bivariate exponential the corresponding point was t/a = .69315. 
That is, the ratio of t, the point at which the effect on reliability 


due to correlation is a maximum, to a, the mean life of a component, 


is constant and is approximately equal to ./7. 


This relation, (4.8), was programmed using FORTRAN language and 


the CDC 1604 computer, and is tabled in Tables 4.01 through 4.15. 


It can be seen that AR(k) is monotone increasing in Pf, and 
hence the value of the reliability will be increased when Pp is 
positive and will be decreased when PP is negative. All tables are 
computed for only positive values of P but because of the symmetry, 
the values are also good for negative P; that is 


LAR(K, - p) = AR(K., p) I. 


There is a separate table for each combination of values of Py 
and Po » for values .995(.001).999. Tables were made for values of 
k from 25(25)1000 and for Pp of .05(.05).25. The ratio of Py to P, 
is designated by s and is shown for each table. Only values of 
Py s Py are used since for any particular case the component having 


the smaller value of p may be designated 1. 


2h 72 





Using the tabulated values, curves were plotted (utilizing the 
CDC 1604 computer) for P, = .995 and Py = .995(.001).999, each graph 
depicts the five values of the correlation coefficient, Pp. [It can 
be seen that the magnitude of the effect increases as the values of 
the parameters approach a single value and reaches a maximum, ee 
P, = Po of .062494 for Pp = .25. Notice that this maximum value is 
approximately .25 ~, which is p times the maximum effect due to 


correlation. (Maximum effect given by Lloyd and Lipow [1] ys 


There is an interesting association between the value of k at 
which the maximum occurs and the "mean'' of the system. If we define 


a System mean as; 


P) + Po 
te | 2 


+ 
ee |= fp) 
2 


* 
and if k is the point at which the effect is maximum, then the ratio 








(4.10) 








de 
k /m is constant and has the approximate value .7. For the case 
where Py = P» the ratio has already been shown to hold (note that for 
that case the mean just defined does in fact reduce to the mean of 


a component). 


Example 4.1 
AS another example, if we take P, = 997. Py = -999, then m = 499. 


Now from Table 4.03 it can be seen that for all values of p the 


function is a maximum at k = 350, and therefore k/m is very close to ./7. 


BAG 





Example 4.2 
If we take the case where P, = coo, P> = .998, then k = 284 


and from Table 4.06 AR(k) is maximum when k = 200, this gives 


k/m = .704. 


This was done for all values of Py and Po which are tabled and 


the value of k/m was in each case = .7 + .O1. 


This relationship might be useful for design purposes in deter- 
mining the component parameter values which would best utilize an 
enhancing interaction between two components, or conversely, spec~ 
ifying parameter values away from these if the interaction is 


degrading to the reliability. 


4.3 Confidence Interval 

To determine a confidence interval for the reliability, the 
independent model was used to first obtain a confidence interval when 
there is no interaction. The method used is that given in Lloyd 
and Lipow ig) page 226. Although other methods are available, see 
Buehler L11], Steck L13], and Madansky [14], there seems to be no 
generally accepted best method and therefore the procedure used was 


picked because of the ease with which it could be applied. 


The procedure for a two component system is to compute 
6 z Ny = fs N, a oS 


N} Ny (4.11) 


and the quantity Na (l - B) = F, 
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Where N. is the number of trials of the fee component and f. is the 
number of failures of the same component, Ni is the minimum N,- The 
number, F, is then considered to be the number of system failures in 
Ni trials of the system. With these as arguments the graphs given 

in [1] page 498 - 502 are utilized to obtain a lower confidence limit 
for any chosen confidence coefficient y. 

Example 4.3 


If we take N, = N, = 1000, y= .95, £, = £, = 1, k = 50, then 
P = (.999)°(.999) = .998001, N= 1000, F= 2. There results a 95% 
lower confidence limit on P of .994. Now since the reliability at 


k = 50 is given by 


R(50) = (P, Py)” (4.12) 


and a lower confidence limit on (P, P>) is given by the lower confi- 
dence limit p; then if R is a lower confidence limit on R, the 95% 


L.C.L. for this example is: 


A $70 


Ro = ¢ ~ (.994)°° = 74 (4.13) 


Bounds on AR(k) for all values of p are * .25 (For proof of this 
statement see [1] page 223) and the limit is attained only when 
p = 7 1. Since p is, however, not known but must be estimated, there 


would also be a confidence interval associated with the estimate. 


A lower confidence limit on R(k) for the dependent model could 


then be given by 


R(k) = R = .25 (4.14) 
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For Example 4.3 the 95% L.C.L. on R(k) would then be R(50) = .74-.25 = 49 
If, however, it were definitely known that P was positive, then the 
Ce. \, R(k) would be .74 since in that case any effect due to correl- 


ation would be enhancing. 


It is recognized that such an L.C.L. on the reliability would 
not be ''good", in the sense of shortest interval, and that perhaps 


a much better technique could be found. 


4.4 Approximating the Effect 
Looking at the relationship obtained for 4QR(K,) given by equation 


(4.8) and examining it part by part, we see that: 


pG +p)G + po) 


VP1 Po 


for values of P) and Po close to l. 





=49 (4.15) 


Taking the next part of equation (4.8) 


k 2k k k 
P -p ap le pe) (4.16) 


It can be seen that this is the product of the reliability 
and the unreliability. Now since the third part of (4.8) is of the 
Same form as the second part, an approximation for the difference 


function can be given by: 


k k k k 
AR (k ) =4o Py o (1 = Py 0) Py 0 (1 - Po 0) (4.17) 
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Notice that this is 4 p times the product of the reliability and the 
unreliability of each component. 
Example 4.4 

As an example of the use of this approximate method, let Py = 
Po * .999 and p = .25, ky = 100. Then using the approximation 
equation (4.17), the reliability difference is : AR (kK) = .007523. 
Now using Table 4.05, the value given for the same conditions is 
2007421. 
Example 4.5 

As another example where this time Py 4 Po> let P, = 2999; 
Py * .998, ko = 200 and p = .10. For this case the approximation 
yields AR (k) = 020568. Using Table 4.06, a value of AR (kK. ) = 


.020543 is obtained. 


It can thus be seen that as a fast approximation the method 


suggested by equation (4.17) does yield good results. 


4.5 Summary 

In this section a jointly discrete distribution, a bivariate 
geometric, was developed from a large general class of distributions. 
It was shown that the particular distribution does in fact meet the 
requirements for a probability distribution function, however, no 
claim is made to its being unique. Using this distribution the 
effect of correlation was examined for a limited range of values of 
the product moment correlation coefficient. Specifically the dif-~ 


ference between the reliabilities when using an estimated correlation 
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coefficient and when the independent model is assumed, was examined. 
Values of the reliability difference were computed and are tabulated. 
A graphical comparison is presented and indicates (1) that the max- 
imum reliability difference occurs when the parameters are equal, (2) 
that the value of the maximum when P = .25 is .25(.25), which is in 
agreement with the known maximum of .25 for P = 1, (3) that the ratio 
a is essentially constant and equal to ./7, ee being the point at 


which the maximum effect occurs and m the mean life of the system. 


A good approximation to the reliability difference can be 
obtained by the approximate method, ie, four times the estimate of 
the correlation coefficient times the product of the reliability 


and the unreliability of each component. 


A procedure for obtaining a lower confidence limit for the 
reliability of a two component serial system is presented, however 
its usefulness is doubtful since it is in no way optimal. No 
better technique could be found although it is believed continued 


research in this area could be very fruitful. 
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COMPOSITE FUNCTION 


5.1 Derivation 

As an extension of the two cases just considered we shall now 
investigate a system wherein the components that go to make up the 
system are continuous and discrete in nature. A simple example 
could be the operating life of an aircraft jet engine wherein the 
operating life depends on the number of starts as well as the total 
running time it accumulates. We may consider the starting function 
as a discrete random variable, specifically the geometric distri-~ 
bution, and the running time as a continuous random variable, 


specifically the exponential distribution. 


With marginal density functions £, (x) = ~ exp (-x/a) and 
f(y) = p” (l-p) we again utilize the theory previously employed to 


generate a bivariate density function of the form 


Egy ov ty ty o)| Me oF) 2) (2Fyl)-Fy(9)-2) (502) 


we sces ce y= UO; elt 2535 


With the familiar restriction that -1 Sv S11 the function fy (x>y) 
is shown in Appendix A.3 to satisfy all the requirements of a joint 
density function, ie, its sum over the range is unity, it is non- 

negative and its marginals are indeed the original density functions 


that went to make it up. The correlation coefficient is evaluated 
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and found to be 


p= vVe/20 + ) (5.2) 


where -l1 Sv S1 andOSp Sl. 


It is again evident that the correlation is restricted to the 
range -.25 $0 S$ .25. In the applications we shall consider, the 
value of p will generally be very high, between .990 and 1.000. For 
this range we see that the value of P is very nearly equal to v/4. 
This implies that for highly reliable items the correlation is 
completely specified by the constant v. Hence we shall consider the 
quantity v = 4 P for most computations. Tables 5.01 through 5.10 


were computed using exact values. 
5.2 System reliability 


The system reliability is a function of the component reliabil- 


ities and can be expressed as 


R(t,k) = P[ x=t, Y=k] (5,3) 


The reliability function is evaluated in Appendix A.3 and the resul- 


tant expression for the system reliability is 


R(t,) = pl exp(-t/a) [ 14v(1-exp (t/a) (1-p*) | (5.4) 
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The system reliability is seen to reduce to the product of the 
component reliabilities in the independent case wherein v = OQ, as 


was to be expected. 


To establish a quantitative measure of the effect of correlation 
on the system reliability a reliability difference function was 
defined as the difference between the system reliability when P = 0 
and that when p #0. This function is denoted by AR(t,k) and is 


expressed as 


QR(t,k) = v pX exp(~t/a) | 1-exp(-t/a)] (1-p*) (5.5) 


The reliability difference is seen to be a linear function of the 
correlation. This function has been extensively tabled in terms of 
the ratio of the total life to the mean life of the components. These 
are denoted by t/a and k/m in Tables 5.01 to 5.10. Further, the 
difference function is plotted in Figs. 5.1 and 5.2 and is seen to 
vary with k/m and is a maximum of .0625 at k/m ®& .693 for P = .25 

and t/a =.7. This is in excellent agreement with previous results 


for the exponential and geometric cases. 


5.3 Confidence Limits 

The subject of deriving confidence limits for the reliability 
function defined above was considered beyond the scope of this 
thesis and was omitted. However, it is evident that this problem 
is of great importance and could well be the subject of a separate 


investigation. 
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5.4 Approximating the Effect 

It is interesting to note that in the event that p is close to 
one the reliability difference function can be approximated very 
closely by the product of the component reliabilities and their 


unreliabilities times 4 ~P. That is, for p™l1 


AR(t 9k) = 46 exp(-t/a) [1-exp(~t/a)] pX (1-p*), (5.6) 


5.5 Summary 

The bivariate density function derived from geometric and 
exponential marginals gives results consistent with those found 
in the previous two sections. Significantly, maximum effects of 
correlation on system reliability occurred for k/m # .7 in all 
cases and the extremum occurred when both k/m and t/a were about .7. 


Again this maximum effect was AR(t,m) = .0625 for Pp = .25. 
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APPENDIX AL 
MATHEMATICAL DEVELOP! (ENT 


BIVARIATE EXPONENTIAL DISTRIBUTION 
To shov fyy(xsy) is a density functions 


must shows (1) fyy(x,y)=0 


(2) , - fyy(xy)dx dy = 1 


-@ -o 


(3) | fy(xsy)ax = fy(y) end 


oo 


fyy (x,y dy = yx) 


—oCO 


where 


fyy(sy) = ~$5 en bVaty/ h + v(1-267*/ (nae) 


(Tt) with aehOg bed} =l=v=1, x20, and y=0 


HG can be seen that <=] = 2e7X/2-S] 
since o=207%/ 2 
therefore since Lave 
i= y(l emer”) 2e 97 = 


and 1 + v(1 we7*/2) (4 ~ 2e°/b) =o 
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now = ge bx/a + y/b)=o 
ab 


since 0 x=e7(x/a + y/b)=1i ana p>0 


Therefore the product 


fyy(ssy) = 4. oW/a + Nha + e(anner*/8)(1-26°9/2)| Bo 


>} [tay ay = | fap * veya + v(1 = 267%/8) 
aa -S — (1 = 20 9/> os dy 


. & 
. mes 


a — [om Paxforv/bay eer [free -20 8A) fo1/P(a-20°9/” ay] 


«oO @ 
look at = | e7X/2qx > ~ | any! bay 
oO Oo 


~ 


Z forx/ax « = Cen =] 





a 


eo 


similarly _1_ i} ebay = 1 
D 


oO 


look at | e7X/8(y ~ 267*/ Vax 
a 


Hep a = C wie X/*) then du = Gow eX/a 


and “for-x- = 0. = 1 for XO , u = =l 
ae he a A 
= = 2 
hence aa fe %/ 8.406 X/ 8a = = | 20 Ga- = =| = 0 
a 
° om] as} 


29 Ja 





Then  E(xy)= [ [* fyy (xy )dx ay 
Oo oO 


Ei) =e gy “ax | Y ov / Davey X grk/8(q26"%/8)ax 
Jy 7 ; 5 


oO 


[z o 3/4 20"9/? Jay 


© 


look at | = 7 X/ “dx: integrating by parts 


o 


“X/Ba~ y= aa eX? 


where u= x/a du = = ax dv =e 


a 


mH 
then | = e7X/8ax = (=) (~a o/8)] = | -= en */Aax 
oO 


° O° 


ao 


[we e*/Pax =O taza similarly [xr e¥/Pay =b 


o o 
Also, if we expand the second term we get 


co 


| x/a e7X/8gx = | dxfa en 2X/agx 


oO 
and if we substitute 2/a for l/a in the interretion by parts, we 
obtain | 2/3. go “ax = a/2 


oO 


Henee from 


E(xy) = | x/a en */ Fa. [ y/b oY Payty i: x/a eM 8/E ay - } 278 e7X/ “a 
CO @ 
| fer aay =x far 2 
O oO 


=—Qheco 





similarly 


a 7 2 
Lf a2 yay= + f[ 2ea= El] =o 
3 =] ool 


Therefore 


| fac dx dy =1el + v[ 0-0] = 1 
o °9O 


(3) 


@ 


/ fyy(x,y)dx = J > gn x/aty/d) 4, 35 oe (x/aty/b) (1267%/A) 
J 4 : (1-2679/ Jax 


oo 


ce orl? | or/2as + ov o¥/P(a~20°0/b) | or/a(r-26°8/8)ax 


“ab 


i 


4, e7¥/P(a) + Ee oA ¥/P(1~2679/b) (0) 


a 
hence [ fxxbsy ee = ot e~y/b = fy(y) 


Oo 


@D 
ee 
similarly | fyy(x,y) dy = "g7 e x/8 = fy(x) 


oO 


To evaluate the correlation coefficients 


We knows E(x) = a, E({y) = b, 0, =a, = 


=95 





we obtain 
E(xy) = ab + v (— “2 ) = ab (1 + w/4) 


Therefore 


_ ab (1 + v/4) = arb 
Oe es 


P= w/t 
Since vies then LOl= ee 
To derive the reliability functions 


R(t) 


P[x=t, yet] I- [ff fypyixey)dx dy 


| Le ax fA +—< oT Payty [+ e7*/8(1-267*/2 Jay 
t t t 


R(t) 


[Re a-20 ray) 
t 


. <O 
look at [+ Se x/ 2... = | | = pate 
t b 


using similar techniques as for the previous integrations we obtain 


Rit) = (e7t/2)(ert/) 4y [ grt/a_gn2t/a)(get/su¢ sae 
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R(t) = eo W/a4t/b) [1 + y (1~07*/8) (1-ert/)] 
To derive the reliability difference functions 
R, (ae qo (t/att/b) [2 nm v(1-e™*/2) (1-37#/b)) 
R, (t) =o (tA + t/d) 
hence since AR(t) = Ry (t) = Ry (t) then 
AR(t) = v e~(t/2 + t/b) (. ~ ort/ay(a — Qt/b) 
“ R(t) = 4 po 't/a #20) — Gat?) weet) >) 
ee eS 
To find the point where OR(t) is a maximum with respect to t/az 


Jet a = b so that 


AR(t) = Ape **/* (1 = o7t/a)2 


Sad =4 Al-2 eH t/2 4 ot / a) 4 ae7et/ 2(je~t/ a)ent/ 20 
a 
hence = (1-e7t/2) + ent/a = 0, ewt/a = oa 

t/a =-ln = 69315 


Max 


9'7= 





APPEIIDIX A.2 
MATHEMATICAL DEVELOPMENT 


BIVARIATE GEOMETRIC DISTRIBUTION 


k 
k+t 
To show that yy pad =l- a 1 
x=U 
k k oo oo 
x X * 
Mepq=a) p =0| Dp = 3p 
x=0 x=O x=0 x=ka] 
Sax ab ; 
note 8 y ———— geometric series and converges 
x=0 ip 
since Pia 
therefore 


oo 


—- x 
he 
x=0 xelekL 
ele hee] 1 
= 5 ae (» + pte - pi*3 + +] 


“te ao po (1 +p +p? + p? “| 


ae 
=e Sop LP | 


aa i=0 


II 
OS 
oe eee 
- 
5 
O 7 
n 
4 
ps 
Wns 





a OL 
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Lemma 1] 
a 
. ; i : 1 ae 
Since the geometric series 2, x” converges to 7 =s~5 provided 


lxi<1, the product a , = can be shown, by multiplying the 





series exoansions, to be 


Se eat oie a tate 
(1x) 


if we miltiply both sides by x we obtain 








= D =x + 2x% + 3x3 + dat + oor 
(1~x) 
Oo 
but this is » ioc 
i=0 
ee ° 
and therefore the series » ix? converges to = D9 provided 
i=0 (1=x) 
Ix] <2 





To show that fyyix,y) is a probability mass function 


when 


fyy (say )=£,(x) f(y) l+y [ 2F (1) -£y(x)-1] [2% y(v)-ty(v)-a 


mm 





where fy (x) = pa Fy (x) ee Gr 


i 
f(y) = 374 oes pv 
with O=p<1 q=l-p OS a 
oe Oe es Se ace VSO lhe 705 eres 


we must show 
(1) fyy (xy) = 0 


(2) YY. tyy(x,y) = 1 
x=O0 y=0 


(3) Pm fyy(xsy) = fy Cy) 9 2 Fn oy) im fy (x) 


x + 
- : | 3 SCR x 
= py "dP 9" dp [a4 [2-2p, “Py "dy taped aa]| 
= x a Sol x yr x 
Py Pa Is a [ epg Py ] [Pa Po” i 


lookine at i= oo a 


=100— 





since 


epee SV aa eee 


x 
= pilot Ane 


oO 
A 


Ll=r(pt1)-122-1 


cS 
t 


=o pee) Sa 


jt 
IV 


then calling 
sor) 
[2 - py “P1* | [2 - pyr - pW] = N,N, 
it can be secn that 


L2=ly h=-1 


Therefore 

ltvly W=0 since -l=v=il_é also 
now since ne qy pe” do =e LOR asc eus Pi» Po then the 
product 


cy aa 
P,"a,P." 9, (1 +v ll I) = 0 





hence 


fyy(xy) = py* a) PoYay (2 + v Ty My) =O 


(2) 


To show that the bivariate geometric distribution sums to 1 


fyy(xsy) = fy(x) £y(y) ty | 28 (x ()-1] [27 7) ty] 


<tr] 


where fy (x) pq Py (x) l-p 


fy(y) = p%q Fy(y) =1-=p 


0 00 a a 
x 
YY Vialw=y > P) 4; P5" 4 Lv [2(1-p)")~py*a,~1] 
x=Q y=0 x=O y=0 


[2(1-p,J ¥1)-p,4 aol] 


f 


a 
# vy x sce] +1 
Py" 9)Pp"dotv|p,%q, [I~ 25 ~P1"a,] p2Ma 1=2p.% ~pJa,] 
2, Ef Pr earl egteloy “| i 1%] P2Vap[ 2 ee 


cw a oo 
= x J ar 2 aCe: 2 < OD 
=). Pa Oy » Pa a [Pp q,~2p G5 -P, ~~ 79 
x0 y=0 eis x=0 eked pe 


2 Ove] 2 
deere eae oe . (A) 
= Oe ape ip 2 ae? 


looking at this term by term 


a 


& 
dpa : ad P oe a = a= 
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looking at 


~ 3 ¢ 2 

y [pa - 2 Gen 0 | 

x=0 

= y oq 2p exe 
x=0 x=0 x= 


a ao 
P< 
1 = 2pq d, (p*)* = o2 Dd (p®)® 
0 x=0 


- at ae ie 

= 1 = 2paq ae since p ab 
i=p 1=p? 

re re | 





Therefore putting these values into equation (A) gives 


Y 


dD fyy(y) = lel tv 00 =1 
=0 y=-0 


fl 


oO 
(3) To show that Py 9) = f(y) 
= 


from equation (A) 
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= r x ot ee 2 
2 far bey) = a py “ay Po” 5 vi [Pa dy 2p] ett et) 
= 


Yq eyrh, me “| 
[Po 9,72Py"  Ag"Py Ap 


heave already shown 


© [pe ae pq" | = 0 


simply substituting these into above equation 


ao 
d, fyy(syy) = py" ay = fy(y) 
x=0 


It ean be seen that by a similar argument 


2 fyy(xy) = f(x) 


To compute E(x,y) and O 


E(x,y) = y y x ¥ fyy(xsy) 
x=0 y= 


from equation (A) and the ecuation preceeding it 


LOL= 





—_— 


au 


y=0 


looking at this term by term 


co oo 
Nixpoa=q | xp 
x=0 x=0 


- wp = pee 
a” a 
looking at 
oo 
Ix] 
> x [pra = 2p q ~ p_2I 
ea) 
2st ax 
= 2 phy 2) xp foe ep ae 


making use of Lemma 1 and the fact p*=1 


we have 


1 Q5= 


= ¥ x 2 as yx 9372p," t04-p 
FYdy 2 J P24 Hn [Py ee niga 


using Lemma 1 
G To ( s ) 


%q47] 
“1 


CO 
Ly [P. Go7KPy InP “9 (B) 





Sy 
2 2 
p , D 5 p 
=— -2q ———- eg OO > 
q (1—p*)* (1-p7)* 
p 2q p? q* p? 
7 q (1<p)*(1+p)* (1—p)"(1+p)* 
2p? pe 
~ g g(ttp)2 (4p)? 
p(ltp)* = 2p2 = qp2 p + 2p2 + p3u2pFap24p3 
q(14p)* q(1t+p)* 
pte? pip) p 
g(1tp)* a(itp)* q(1+p) 


Putting these values into equation (B) givess 








P 
won (2][4 
o(14p,) a, (14) 


Py Po Vv 


ee 
Gy Ay (1+p, )(1+p,,) 


To conpute Q 





co qp p 
E(x) = Lx we i ep ea 
x=0 (1p) q 


Making use of Lemna 1 
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p 


. Zo a ee eee : 
Lioyd & Lipow /2| p 123 gives C,, = for this distribution, 


a” 


p 
therefore on =) 


now using the defining equation for Pp ae, 


E(x,y) — E(x) Ey) 











Ox Oy 
elves 
1+ ——|-—- — 
G4 Io (14p, )(14p,) d, qo 
p = 
V Pr V5 
oa Lo 
P2 Pe eae ae : 
i oe a 
* "2 [(a4py) (4p) 
p= P1 £2 SaNice” ave = ok 





ra (1+p,) Hoda, <4 


then ~ z= (SEs 
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To compute R(k,) 


Rk =P [X2kos ¥ =k] = za yf f(x,y) 
— om 


‘. yo . ey 2) ty(y) (240 [ 2F (ac)~t.(x)=1] 
x=k Ban 


2, or-£yy)-1]]] 


now f(x) = 5g : F(x) = 1 ek 
fyly) = pq : Fy) =1- prt 
Therefore 
R(k, 3° y EB GPa do bt [2(1p, ) pay os 1| 
xk, Y=k, 


2p") -p Jay i} 


looking at this term by term 


oes 





As 
~ 
13 
m4 
52 
(scat 
Oo — 
6 
ore 
© 
+ 
- 
pes 
ms 
oO 
ok 
nN? 
o 
@ 
gs 
xno 





Pq-@ ee 
xk, AK 
k : 
= a) Pee eee tp? + evel 
Keo Ko i k. 
ee ie chi 7 ae 
L=0 L-p 
looking at 
— : J ca) « 
s [ aptq 2p hg yp 9g” “pa | 
x=k, 
oo oO 
x e Poth Pp ane) 
=) pa- 2), D a=) p*q 
2=Ko X=Ko xek 


i ee a ee 
p ~ 2pq ), (p*)* = ° a eee 
Do, 
© 


7 
som 











2 k ak 
k (Gre ; (p~)yo 
a. i. <pq ? ee ry 
leap” Jens 
oD a: 
by putting (p<) Loo Oe >. a = 
mae Lo 
Kk Ie 
ee Oe) dee 
itp 1+p 


| Nef] a 1, 
p(1tp) ~ 2p Kot . phe + peKotl, 


iI 5 
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Ae aD 


(1+p) pm pro (1+p) 


1l+p 

= pro eo =o - ene) 
Therefore putting these values back into equation for R(ky), 
and using the corresponding terms for y, gives 
Rs.) = py Py PH | [ee Py?) ] [ Pz" = P,) } 
iene Call tiie Ro (ky) and call R (ko) the reliability obtained 
when using the product rule, is. 

Ry (kg) = py? pp 


then the difference, which we shall call AR(k,) becomes 


R(k,) 


Ry(ky) = Ry (ie) 


Vv [m° (1 = 1°) | [ p,*° (1 = Po?) 
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now putting it in terms of Pp gives 
(1+p, ) (1+p5) 


VP] Po 





1 . 
AR(k,) = p lp, ‘0 ( 1-p,“° ) p> (1=p se ) 


To determine the value of ky at which AR(k,) is maximum 
a Ko Ko ko “Ko | 
OR(k,) =v [Py Py Py P, 
if we take the case Py = Po then 
k 


we seek the derivative of this with respect to ae 


(v is independent of k,) 


d 4R(k,) 





Ky a d 
dk, ak 





[Pp - p*?} 


k 
2D Ono) (porn p- 2p"O1n p) 


k. 2k k 
2 v (p Pp?) In p (p © = 2p**o) 
setting this equal to zero and solving for i 
ic, 2k Ky ak 
2vinp (p 9% p™°)(p ° — ap"9) = 0 


ld 





(po - po) (pk = 2p~Ko) = 0 


one or both of these terms must be O3 if we take the first term 


ro : po =, 
es = pro 
Cee a 
k =2k 


& a6 ke, = 0 
and it can be secn that this is a minimum point. 
Now taking the second term 


oo 2 2 pao 


i 
© 


k 
ee eae 


ee Dl 2 2 ae 


k Gi p= o> np v= ine 


ee 





dene 2 in 2 








ie SOS 
inp=<-2i1inp lp 
hence 
Ee 
Ki fee 
. dbelsye 


is a maximum point of AR(ky) when p, = Po. 
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APPENDIX A.3 
MATHEMATICAL DEVELOPMENT 


JOINT EXPONENTIAL, GEOMETRIC DISTRIBUTION 
To show that fyy(x,y) is a probebility function 
when 
PryOoy) = yx) yy) [ I4+v (2F (e)-1) (aR (y)-1)] 


where 
-x/a 


a 
Paces ac Om O=<ax=co 
4 a 


fy\y) =pq y=0, 1, 2, o0* 5 0=p=l 


with q=lep , -l=v=l 


a 


| £ Cx! Jax’ =1 - ew/a 


Fy (x) 


0 


J i 
FL(y) = \ pPiqzl-p 


This gives 
1 -x/a, Vv ox /2, al yY 
Ley) = —e wqlityv [2(1 <1 oo fea - p!")=p ql (B) 
a 


-l14= 





Since fyy(xsy) is a continuous function over a countably infinite 


number of points, ie. 


Z, 
f f(y) dx has a value at integral values of y only 
xt0 


some positive value if y = C0, dye 2, hese 


Ey %y) qx = 


0 otherwise 
x=0 


Then to evaluate fyy (x,y) over its entire domain we must look at 
an infinite sum of such integrals. 


To show that this fyy (x,y) is a probability function we must show 


(2) x ” ferley) = 1 
x=0 
J 00 
(3) 2 fyy(xsy) = fy(x) , fpylery) dx = £,(y) 
x=0 


oo] 1 5 = 





(1) using equation (B) 


i 
Eyy (xy) — /a pq | ltv Be = o7*/2) = 1] 


[2(1 Se) pas 1] 


Be veeee E = 207 *) E ~ pytt . py 


a 
In Appendix A.l it was shown that 
sigh ee a 
and in Appendix A.2 that 
~1=1 - ee - py =1 
hence the quantity 
1 +v/1 - 20/8 Jf a ~ pytl _ py] =o 


, zi x/a 
now since =e =o 
a 


and py G0 


ale 
the product et i py G0 
a 
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Therefore 


ee ee) oe 
a oy) 





oO ce 
(2) To show » fyy(x,y) dx = 1 
y=0 


x=0 


= _ = r a ~x/a -x/a =x/a 
Xe f Cy bv or 2 fk e platy [e (1=2e ) 
x=0 


x=0 
ate 
pq (1~2p) tap of 
a 
© at 3 co 
=) pq aes eee Sp g(1=2p) tpg) 
y=0 a. y=0 
x=0 aed 
: ool / ‘| 
x=0 
= (1) (1) tv /0*0 | =] 
calculations follow, 
Oo 
us ce 
to shows —e */84.. aL 
a 
Oo 
ce co s 
— das = [ne a dx) == owe |= - [0-2] =o 
a 
Oo ° O 
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to shows = 





orca 
y=0 
co 
>. pq = aS p” geometric series which converges to 
y=0 y-0 
if |pl<1 
Ie es 
Sone i: 
= 0 ee since q=l-p 
=e 
@w 
to shows 
: e 2/84 = ae X/®) a= i0 
O 
© a 
O O 
oD ©0 
-x/a 1 . 5 
=-aj © pane g) [eZ an 
0 O 
0 a 
ae owe ta ee = -2[0-1] +a [ 0-2 ] = a-a = 0 
O O 
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to shows . J ao 
2 Pall = a - pq) =0 


ao 

Se aries ee ae rae 

y=0 LP a~ Epi 
y=0 








since pr<=l 


2pq q* 
(1=p)(1+p) (1=p)(1+p) 














ee 2p i“ q 
LS Lio 
ae —2p = (1 = p) 
foe 
Ot 
=j]+ ie ee 
neers 
=p <—1 
=o 
1+p 
- Lote 
= on =] = = 
nde 1 0 





(3) To show J f(xy) = f(x) 
y=0 XY Xx 


co rae) ale ; 
Y fpr) = epg [1 [2-07/*) ~ 2] 
y=0 


y=0 
[2(1—p9"*)=p¥ =] 


on : pea ty [ o/ 2(1-26°3/ *] 
y=0 


e [eu 


[ & eat ROY a pa) 


but we have already shown,in Appendix A.2, that 


and that 


oo oo 
ye a(1 - apg) = ¥ [ P% - peer VQ? | = 
ee 


hence 


{1 
| o 
@ 
ve 
Be. 
© 
vV——_ 
kK 
+} 
< 
— 
o, 
a 
a) 
ro 
q 
~w 
MO 
4 
os 
09 
=, 
—, 
Cy) 
ag 


: 
fy) 
ay etl 


| 
° | 
es 


~120= 





co 


To show [ fy (x,y )dx = f(y) 
x=0 


o0 = in : 
/ Pyy (x,y )dx = — pa ey [2(a-6/*)-1] 
x==0 ae 
[2(1-" +1) py q-2] dx 
= fe [ae pene ee en 
x=0 x=0 


but we have already shown, in Appendix A.l, that 


a0 
1 
—X 
i [2 ... =] 
a 
x=0 
Oo 
ene that == = ne7*/2) ax —6 
x=Q 


hence 


aD 


[favor = p’q(1)4v p%q [2(1~p" T=) a2 | [0] 


x=0 


y 
pq iV) 
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To calculate EK,y ) end Q 
Here we must again look at an infinite sum of integrals of 
Ey (sy) weighted by the function xy. 


CO (e @) 
E(x,y) = ys [» fyy (x, y)dzx 
y=0 


x=0 


from equation (B) 


= i 1 =x /a 
E(x,y) = ) J x —e / yp | liv [2(a-e*/*)-2) 
y=0 a 
x=0 


[2 (1=p¥"+)-p¥ eI] ax 


ao a a 

Dy Oinig 3 yp’o [= e7X/ Bg a ats i= x @ 2/8¢9 997%/2)a5 

y=0 a a 
x=0 O 


y yp a(1-2p" ~ pq) 
y=0 
co 
looking at | a ee 
J a. 
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using integration by parts: 


[rev uv - [ve 


x mS 
Tete . seen ae 
a 
a 
then du= —— dx ve fear =~ 
a 
oO ao co 
ae ee lt) 
gives: e X/a So vs] or. ae 
a a e a 
O O 
w 00 
wr / x Be 
= — xe %* w26 me == ¢ */8(ta)] 
O O 
ni B 
a) = tin -<9/*(ta)) Sr [-<78/*(a4a)40°/ *(ota)} 
O ae co O B—>o0 
= lim fe - oB/ “(Bta)] 


B-,oa 





= ne eas dm, B Rh eo lim a oo 
B-Sc° Bo}@ B—»,co 
B 
=a-a- lim =- 0 
B-o ob/e8 
B 
=a- lin 
ois 9/2 
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B 
using L' Hopital's rules lin Ba = Jim 
e 


B00 B 00 d _b/a 
dB 
ie 
= lim : = 0 
ee =b/a 
Therefore 
co 
ae 
mX/ a. 
—e ay dx =a 
a 
O 
@m 
x 
js om”: a @ a 
looking at —Bes fa (1 = 26 / ydx 
a 
° 
aD ce 
x ex 2x/ 
= ——6¢ / ax —-j{—<e dx 
a a 
o 
first integral is a (from above) 
Z u 
second intesrels if we substitute — for —=- in above, we get 
a. a, 
@ 
OX a 
=2x/a. 
= 6 ax = —— 
a a 
O 
ao re @) 
x am a a 
: aX / 2 2X / a 
a a. Z a 
O 


~l2h= 





cD 
looking at >. y p q 
ae 


malcing use of Lemca 1, Appendix 2, ve see that 


go p 
\ yp’ = 5 since Ipl< 1 
y=0 (1-p) 





hence 








: -y ° o Pp pg p Pp 
yrqtaq ypraq ee Tay ose 
y=0 y=0 (Q-p)*  (2-p)(1-p) ssp q 
= 2ytl 
now look at >. (~2 y p q) 
y=0 
@o p* —e pa 
= 9 S ( Zs As Se we Le 
y=0 (1 = p*) (1=p*)(1-p<) 
-2p?q - 2p? 
(lep)(1itp) (1p) (1+p) q(1 + Be 
now look at 
2 Die 
00 oO p =~ q~ p 
ey pq") = ~o~ > y(p")¥ = ae ee ae 
y=oO y=0 (1=p*) (1=p)(1+p) (1~p) (1+p) 
= e 
(cto) 


| 250 





co 
Therefore y pq (1 — 2pvtt pq) 
y=0 


ao 


= Ly wa- aby a= by wae 


y=0 
= 2p? p* 
q ality (py 


p(1+p)* = 2p? = gp* 
q(1+p)* 


p(1+2ptp2) - 2p? - (1+p)p- 
q(1+p)* 


pt2p? + p? = 2p? = p* + p? 





q(l+p)* 
7 p*tp - p(1+p) ee: 
q(14p)* q(1tp)* q(1tp) 


Therefore putting these values into equation for E(x, Y) gives 


p 
ne | | 
al 
an S Vv 
1 + 
q 2(1+p) 


a 


2 


Pp 
(a) 





E(x, Y) 

















q 











To find expression for PD 














p “fp 
We have EG. » O, =a, af a aes, 
q q 
E(x,y) ~- E(x) Ely) 
and train 
Ox Oy 
hence 
ap V ap 
[2 . | - . 
, Ne 
q 
V 
: (2 * 2@4p) 7 1 
Vp- 
v Vp 
i) cae lvisel O<=p<l 
2(1+p) 


therefore | fe) =} 


To evaluate R(t,k) and determine AR(t,k) 


R(t,k) =o | X=t, Ek | = >: , £ (x97) 
yr-k 
x=t 


where again we must sum over y and integrate over x. 


~l27= 





oD 
1 

R(tyk) = [— on Fay Ded pb. g. Tr {fe ¥/2(4 ~267%/2 ax 
a 


+1 
5 p¥q (1-2p7""=p"q) 


y7k 
looking at a oO 
1 ux/a as = 
rere m4 ie 
a t 
t 
looking at = 0 go 
au cy 2 
— 97%/8(1~26 1/8) a = |—e */ A. ~{—e 22/84 
a a, a 
ie 
2 aL 
from above,first integral is g-t/ . substituting —-for —éin above 
a a 
co 
2 
integral gives fi ax = e72t/a 
a, 
@ 
therefore |; 1 
fi #/8 05 ng -2/2 Jax = o7t/ ae t/a = ~t/ 2(y-67t/ = 
a, 


t 
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looking at 


+ pktl + pt? + +e 


=qp* [1 +p + p> + pF + ere 


aE 


I 
nn 





aP 
1-p 


looking at 


CO 
», Pa(2 — 2p" = p¥q) 


y=k 
oo Oo cD 
yr-k y=k yok 


cO @ 
po -2 pq Yi (p27 =o D (p*)¥ 





y-k yk 
k k 
(p*) (p*) 
meb = <bG oo 
1 eae 1 = p* 
2k+1 
_ ok ap (1 - p) p™* 
Lap 1 tp 
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an sales oprctl ee 
Iva 
k + k+l 2k 
p + p fe Pp =p p“(14p) ~ p“(14p) 
1a (1+p) 
k 
= p* ~ p= p* (1 = ph) 





Therefore putting these values into equation for R(t,k) cives 
R(t,k) = eta pk + y er*/a(a ~ en t/a) + pk (1 = P| 
calling this Ro(t,k) and defining R (tsk) as 


R,(t,k) = e7 t/a pX (product rule) 


Then if AR(t,k) Ro(t,k) = R, (t,k) 


ARGC) 


vo t/8 ok (y wo 8/2) (q — pk) 


or in terms of 2 
2 P(1+p) 
Yp 


Ano) e~t/apk (Gre ae p') 
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